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In this part of the paper, the higher order approximations to the general dispersion 
equation are considered. It is shown first that the results of the first approximation are 
also found to be valid in the neighborhood of Q = (), even though the phase velocity goes to 
zero for that value of $2. The second approximation which gives correctly the phase veloc- 
ities of the order of the acoustic velocity a, in the electr on g as or lower and the third ap- 
proximation which yields the phase velocities of the order -y/aCo are worked out. A synthesis 
of the results of the first three approximations yields the dispersion curves for the three 
independent modes which propagate in an electron plasma. 

The dispersion curves for the three modes obtained as a result of the first three ap- 
proximations are not valid in the axial boundary layer, where the dispersion is found to 
undergo rapid changes in the close nieghborhood of the frequencies il= R and 12=1. Using 
the boundary layer techniques, the dispersion equations in the axial boundary layer, for 
the frequencies in the close neighborhood of 12= R and 12=1 are obtained and analyzed. 
The boundary layer analysis clarifies the mechanism of coupling between the purely longi- 
tudinal plasma waves and the purely transverse electromagnetic waves. With the knowl- 
edge of the above coupling that takes place in the axial boundary layer, the three inde- 
pendent modes may be properly designated using a procedure which is found to be uniformly 
valid for all directions of propagation and the inconsistency found in the magneto-ionic 
approximation is removed. 



1. Introduction 

The characteristics of a plane wave propagating 
in an unbounded, collision-free, compressible elec- 
tron plasma along a direction which makes an 
arbitrary angle with that of the static magnetic 
field are studied in this paper. In part I of this 
paper, the cubic equation which specifies the square 
of the phase velocity of the plane wave was de- 
rived and was shown to factor out into a linear 
and a quadratic equation for the two special cases 
of propagation, namely those along and across that of 
the external magnetic field. This feature enabled 
the dispersion relations for the above mentioned 
special cases to be obtained exactly without any 
difficulty. For the general case of propagation in 
an arbitrary direction with respect to that of the 
static external magnetic field, the exact dispersion 
equation lent itself to be approximated by a quad- 
ratic equation in the situation where only the 
phase velocities of the order of the velocity C 
of the electromagnetic waves in free space or higher 
are sought. The analysis of the dispersion equa- 
tion obtained in the first approximation yielded 
zero values for t lie phase velocity for the three 
frequencies, 12=0, 12 3 and 12 4 . Consequently the 
validity of the first approximation becomes ques- 
tionable in the neighborhood of the above three 
frequencies. 



In (his part of the paper, it is shown that the 
first approximation is accurate in the neighborhood 
of the frequency 12=0 despite t lie fact that the 
phase velocity becomes zero for 1 hut frequency. 
A second approximation to the dispersion equation 
is obtained which gives accurately the phase 
velocities of the order of the acoustic velocity a 
in the electron gas or lower. The second approxi- 
mation yields a linear equation and its analysis 
shows that the phase velocity becomes infinite for the 
same two frequencies 12=12 3 and 12 4 for which the 
first approximation yielded zero values. The first 
two approximations are obviously not accurate 
in the neighborhood of the two frequencies 12=12 3 
and 12 4 . A third approximation which is valid 
in the neighborhood of the two frequencies 12=12 3 
and 12 4 is obtained and it shows that the phase 
velocities for those two frequencies are of the order 
•y/aC d . Moreover, the results of the third ap- 
proximation are found to merge continuously with 
those of the first two approximations, thus yielding 
the oa erall dispersion curves for the three inde- 
pendent modes that propagate in an electron plasma. 

The dispersion curves for the three modes obtained 
as a result of the first three approximations, change 
as the direction of propagation is changed and 
attain the known values for the propagation across 
the static magnetic field but not for the case along it. 
The reason for this behavior is shown to be due 
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to the fact that in a small region about the direction 
of the static magnetic field, called as the axial 
boundary layer, the dispersion in the neighborhood 
of the two frequencies ti = R and 12=1 changes 
very rapidly. By employing boundary layer tech- 
niques the dispersion equations in the axial boundary 
layer in the close neighborhood of the frequencies 
Q = Ii and 12=1 are derived and analyzed. This 
analysis enables the proper understanding of the 
mechanism of coupling between the purely longi- 
tudinal plasma waves and the purely transverse 
electromagnetic waves for the case of propagation 
other than that parallel to the static magnetic field. 
Once the mechanism of coupling is understood, it is 
possible to designate the three independent modes 
that propagate in an electron plasma using a 
procedure which is uniformly valid for all directions 
of propagation and the inconsistency found in the 
magneto-ionic approximation is no longer found to 
exist. 

2. Higher Order Approximations for the Re- 
gion Exterior to the Axial Boundary Layer 

It will be shown now that even though the phase 
velocity becomes zero for 12=0, the approximation 
(25), which was derived from (24a) on the condition 
that the phase velocity is of the order C or higher, is 
still valid in the neighborhood of 12=0. In the limit 
of 12 tending to zero, (25) together with (24b, c) 
become 



where 



\*+Ajp\*+AP=0 (30a) 1 

2 1 (1) = 12 2 Q2+A >2 (1 + / 2 )] (30b) 

2i 1) — tfCffiH 2 . (30c) 



If 12 is of the order a 2 1 CI, the solution of (30a), on 
omitting terms of the order a 2 Id in comparison with 
unity, is obtained as 



xW-l^ar 2 /?/. 



(31) 



For 12 of the order a 2 1 CI, (31) shows that the phase 
velocity is of the order a. In the limit of 12 tending to 
zero, the exact dispersion equation, after omitting 
the terms of the order a 2 1 CI in comparison with unity, 
becomes 



\*+A?>\*+AP\ 2 +AP =0 



(32a) 



where A\ {1) an d ^ 2 (1) are the same as in (30b) and 
(30c) respectively and 



A^ = -tfCta 2 RH 2 . 



(32b) 



When the phase velocity X is of the order a and 12 is of 
the order a 2 /C%, on omitting the terms of the order 



d 2 ICl and a A /Co in comparison with unity, it is found 
that only the first and the third terms, which are of 
the order a 6 , need be retained on the left side of (32a) 
with the result 

(33) 



x 2 =V-^>. 



Note that (33) obtained from the exact dispersion 
equation (24a) is the same as (31) which was obtained 
from the approximate dispersion equation (25). 
Therefore, the validity of (25) in the neighborhood of 
12=0 is not impaired even though the phase velocity 
goes to zero for that frequency. It is to be noted with 
the help of (31) that the phase velocity goes to zero 
as 12 1/2 in the limit of 12 tending to zero. 

If the phase velocity is of the order a, it is seen 
from (24a) that the fourth and the sixth terms on the 
left side of (24a) are of the order da 2 and the re- 
maining terms are at least of the order a 2 /C% lower 
than these two terms. Hence, only the fourth and 
the sixth terms need be retained with the following 
result 



\ 2 = 



a 2 12 2 (£2 2 -W) 



AP &-Q 2 (l+R 2 )+R 2 l 2 



q 2 12 2 (12 2 -ffT) 
: (12 2 -12i)(12 2 -12 2 )' 



(34) 



1 The equations and the figures in parts I and II of this paper are numbered 
consecutively for the sake of convenience. 



Since !23<Ci? 2 / 2 <C^i (34) gives real phase velocities 
only for fl§< U 2 <RH 2 and 12 2 4 < 12 2 < oo . In figures 6a 
and b, the dotted portion shows the real phase 
velocity given by (34). In the second approxi- 
mation, the phase velocity diagram is seen to 
have two branches, one of which goes to zero for 
12 2 =i? 2 / 2 and the other asymptotically approaches the 
acoustic velocity a in the electron gas as 12 tends to 
infinity. Since the approximation (34) is valid 
strictly only for the phase velocities of the order of a 
or lower and since the phase velocity given by (34) 
becomes infinite for 12=12 3 and 12=12 4 , it is obvious 
that (34) also, is not a valid approximation in the 
neighborhood of the frequencies 12=12 3 and 12=12 4 . 

A third approximation to the exact dispersion 
equation (24a) can be obtained in the following way 
to yield phase velocities which are reasonably 
accurate for frequencies in the neighborhood of 
12=12 3 and 12=S2 4 . If the phase velocity is of the 
order ^aC , it is found that the terms on the left- 
hand side of (24a) are of the order 

aC K& aCl (S) aCl (&' 

aCl, aCl (jjj> and aCl (0 

respectively. Hence only the second, the fourth, 
and the sixth terms need be retained with the result 
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The solution of (35) is 
A « 



2M 



± 



VG#)~ 



(36) 



From figure 4, it is evident that A ( $<^0 in the 

neighborhood of & 2 =&3 and that ylio^O in the 

neighborhood of 12 2 =124. Also from (24c) and the 
fact that S2r;</« 2 <1^, it is clear that A { &>0 in 

the neighborhood of 12 2 =12l and A$<0 in the 

neighborhood of 12 2 -^. Hence A { £/A^<0 for 
the value of 12 near both 12 3 and 12 4 . Therefore 



\V2i!jv ~ 



^-30 

AW 
^■10 



> 



9 /4 0) 



and only the upper sign in (36) will lead to a real 
value for the phase velocity and, consequently, the 
lower sign may be disregarded. 

When 12 2 is sufficiently removed from either ill or 
QJ, it is found with the help of (24b-d) that A™/ 
Aft is of the order a 2 /6 r2 . less than (A { i ) ) 2 /(2A\\}) 2 . 
The square root in (36) may, therefore, be expanded 
yielding the following result: 



A (] 



2 A !,V 



\\A{ 



(37) 



For 12 2 >12 2 , it is seen from figure 4 that A$<C0 and 
A2K0. Similarly for 12 2 >12 2 4 , A§>0, and 48> } >0. 
For ft 2 sufficiently greater than £2 2 or Bj, A®/ 
A { io ^>0 an d, therefore, (37) becomes 



/1U) 
A- T7TT 



(38) 



It is found, on comparing (34) and (38), that the 
phase velocity (38) obtained from the third approxi- 
mation becomes identical with that given by the 
second approximation for 12 2 sufficiently greater than 
S2| or fij. 

For 12 2 <12 2 , A ( £<0 and A%>0. Also for 12 2 <12*, 
^io>0 and A%<0. Hence for 12 2 <13 2 , and 12 2 <r2 2 , 
on omitting the terms of the order a 2 \C\ in com- 
parison with unity (37) becomes 



X 2 =- 



AW 



(39) 



As 12 approaches from below the close vicinity of 
either 12 3 or 12 4 , it is found with the help of (24c) 
that the second term inside the radical in (26) 
becomes very much smaller in comparison with the 
first term, enabling the radical in (26) to be expanded 
yielding the following result: 



\ 2 - 



2A 



, \A8?\ A® 
r 2A ~ i ~\Att\ 



(40) 



The upper and the lower signs in (40) correspond 
respectively to the extraordinary and the ordinary 



electromagnetic modes. For the O 11 mode and for 
12 slightly less than 12 3 , figure 4 shows that A { £<C0 
and this reduces (40) to 



XS= 



AW 

^20 

" A &)' 



(41a) 



Since, as seen from figure 4, A'^^0 for 12 2 < G|, (41a) 
leads to real values for the phase velocity. In the 
same manner, it is seen from figure 4 that A { ^^>0 
for the e 11 mode and hence for 12 slightly less than 12 4 is 
found to yield 

X^=-4S- (41b) 



Using figure 4, it is found that -4g}<0 for £2 2 <f2 2 
and therefore (41b) gives real values for the phase 
velocity. A comparison of (39) with (41a, b) shows 
that the phase velocities obtained from the first 
approximation for 12 2 very close to but less than either 
fif or 124 are identical with those obtained from the 
third approximation. 

Also the third approximation (35), as seen from 
(36) and the argument following it, gives the follow- 



ing expression for the phase velocity for !2 2 ^ 



-v ; 



MO 



--tii and 



(42) 



Since A^/AftKO for 12 2 =12§ and 12-= 12^, (42) is seen 
to yield real values for the phase velocities for 
12 2 = 12r; and 12 2 =13i In addition, with the help of 
(24b, d), (42) is found to give a phase velocity of the 
order -\aC for 12 2 =12 2 and 12 2 =124. Consequently, 
the third approximation (35) which was derived 
from the exact dispersion equation (24a) on the 
condition that the phase velocity is of the order 

•\aCo, i s certainly verv valid in the close vicinities of 
both Q 2 =Q$ and 12 2 =12t 

The synthesis of the results obtained in the three 
approximations shows that the dispersion curves 
for the () l and the e 1 modes obtained in the first 
approximation and therefore strictly from the 
application of the magneto-ionic theory are practi- 
cally unaltered by the two higher order approxima- 
tions which, in essence, include the effect of the 
compressibility of the medium. However, the n 
and the e 11 modes are affected in the following 
manner. The O u mode as given in the first approxi- 
mation, instead of going to zero for 12=0 3 , continues 
and goes to zero for Q?=RH 2 . The dispersion for 
the 11 mode is given by (29) for 0<12 2 <12^, by (36) 
in the close neighborhood of 12 2 =12?< and by (34) 
for Ul<CU 2 <RH 2 . On the other hand, the e 11 mode 
as obtained in the first approximation is modified 
by the higher order approximations in such a man- 
ner that the phase velocity instead of going to zero 
for 12=12 4 , levels off to the value of the acoustical 
velocity a in the electron gas and thus continues up 
to 12= 00 for which value the phase velocity is 
exactly equal to a. The dispersion for the modified 
e 11 mode is given by (29) for a?<fi 2 <$, by (36) in 
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the close neighborhood of tt 2 =&l and by (34) for 
124<12 2 <oo. The figures 6a, b are respectively re- 
drawn in figures 7a, b wherein the third approxima- 
tion (36) is included for the phase velocities in the 
close neighborhood of Q 2 =tt 2 3 and 2 =0* and thus the 
complete dispersion curves for the arbitrary direction 
of propagation of the plane wave are obtained. 

It is convenient to redesignate these modes in such 
a way as to fit in with the existing nomenclature for 
these modes for the two special directions of propaga- 
tion, namely those along and across the static mag- 
netic field. The 1 mode is obviously the modified 
ordinary electromagnetic mode [MEM(O)] since it is 



10 r 



. MEM(O) 



: '/2 



1/2 




Figure 7a. Complete dispersion curves for the arbitrary 
direction of propagation (R 2 =l/2, Z 2 = l/2). 



cutoff for Q=l, asymptotically attains the velocity 
Co of the electromagnetic waves in free space and 
has always a phase velocity greater than C . The 
modified e 11 mode is redesignated as the modified 
plasma mode (MP) since its phase velocity asymp- 
totically approaches the acoustical velocity in the 
electron gas. The e 1 mode is redesigned as the modi- 
fied extraordinary electromagnetic mode [MEM (J? 1 )] 
since its phase velocity also asymptotically ap- 
proaches Co in the limit of infinite frequency. The 
modified 11 mode for large values of B may overlap 
with either the MP or the MEM(O) modes and 
therefore, should be another branch of the modified 



R 2 = 4/3 



f-\/2 




Figure 7b. Complete dispersion curves for the arbitrary 
direction of propagation (R 2 = 4/3, I 2 = 1/2). 
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extraordinary electromagnetic mode [MEM (X 11 )]. 
The reason for naming the second modified electro- 
magnetic mode as extraordinary is perhaps due to 
the fact that it has a branch wherein the phase 
velocity has considerably lower values than C . 
Also unlike the other modes, the first part of the 
MEM(A MI J mode has a phase velocity that increases 
with frequency. 

In the absence of the external static magnetic 
field, the plasma and the electromagnetic modes are 
uncoupled. The plasma mode is entirely longitudinal 
and the electromagnetic mode is totally transverse. 
However in the presence of the static external mag- 




Figure 8a. Complete dispersion curves for the arbitrari 
direction of propagation (R 2 = l/2, I 2 = 0.95). 



netic field, the modes of isotropic plasma get mixed 
so that all the modes in an anisotropic plasma have 
both longitudinal as well as transverse field compo- 
nents. Hence all the three modes in an anisotropic 
plasma are designated as modified in order to be able 
to distinguish them from the corresponding modes in 
an isotropic plasma. 

Even though the three modes in an anisotropic 
plasma have both longitudinal and transverse field 
components, it is found that each of the modes is 
predominantly longitudinal in certain ranges of fre- 
quency and predominantly transverse in the remain- 
ing ranges of frequency. The modified ordinary 
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Figure 8b. Complete dispersion curves for the arbitrary 
direction of propagation (R 2 = 4/3, Z 2 = 0.95). 
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electromagnetic mode is predominantly transverse 
throughout its propagation range 1<C0<C°°. The 
modified extraordinary electromagnetic mode is pre- 
dominantly transverse in the frequency ranges 
0<CO<0 3 and 2 <0<Coo and is predominantly 
longitudinal in the frequency range 3 <0<Ci^ and 
in the neighborhood of 12= 12 3 , the rapid transition 
from a transverse to a longitudinal wave takes place. 
The modified plasma mode is predominantly trans- 
verse in the frequency range Oi<CO<0 4 , rapidly 
changes character in the neighborhood of the fre- 
quency 0=0 4 to become predominantly longitudinal 
in the remaining frequency range 4 <0<C °°. 

As the direction of propagation of the plane wave 
is changed so as to approach the direction per- 
pendicular to the static magnetic field, I 2 becomes 
smaller and this has effect only on the values of 
3 , Rl and 4 . It is found that I 2 approaches zero 
both 3 and Rl tend to zero and 4 tends to the value 
V1 + -R 2 . It follows therefore that as the propaga- 
tion direction is changed so as to approach the 
direction perpendicular to the static magnetic 
field, the MEM(O) and the MEM(J J ) modes remain 
practically unaltered whereas the propagation range 
of the MEM(X n ) mode gets continually diminished 
until it finally vanishes in the limiting case of 
propagation across the static magnetic field. For 
the MP mode, the only effect is the frequency 4 
around which it changes from being predominantly 
transverse to predominantly longitudinal, continually 
increases as I 2 is decreased and approaches the value 
of -y/l-\-R 2 in the case of propagation across the 
static magnetic field. A comparison of these 
results with those obtained for the special case of 
propagation across the static magnetic field, as 
depicted in figure 3b shows that the dispersion for 
the three modes obtained here for the case of prop- 
agation at an arbitrary angle to the direction of the 
static magnetic field, attains the proper limiting 
values for propagation across the static magnetic 
field. 

On the other hand, as the direction of propagation 
is changed so as to approach that of the static 
magnetic field, I 2 increases continually and attains 
the value unity. As I 2 increases to the value unity, 
Rl approaches R and 3 and 4 approach respectively 
the smaller or the larger of the two values, 1 and R. 
It is found that the dispersion curves for the propa- 
gation direction at an arbitrary angle to the direction 
of the static magnetic field do not in the limiting 
case of l 2 =lj become identical with those for the 
case of propagation along the magnetic field as 
depicted in figure 2b. Consequently the results of 
the higher order approximations obtained in this 
section hold good everywhere except in a small 
boundary region about the direction of the magnetic 
field and, therefore, are said to pertain to the region 
exterior to the axial boundary layer. In the axial 
boundary layer, the dispersion curves undergo 
rapid changes due to the decoupling of the plasma 
mode. A different type of higher order perturba- 
tions is therefore called for in the case of propagation 



very close to the direction of the static magnetic 
field. 

The dispersion curves for the three modes are 
drawn in figures 8a and b for two values of the param- 
eter R, namely (i) R 2 =y 2 and (ii) R 2 =% for the case 
in which the direction of propagation is very close to 
that of the static magnetic field. Specifically the 
value Z 2 =0.95 is used. An examination of figure 8a 
shows that for the case 7? 2 < 1, the dispersion curves 
for the modified plasma mode and the modified 
ordinary electromagnetic mode approach very close 
to each other in the neighborhood of the frequency 
0=1. Similarly from figure 8b, it is seen that for 
the case 7? 2 >1, the dispersion curve for the modified 
plasma mode approaches both the curves for the 
modified extraordinary and the modified ordinary 
electromagnetic modes in the neighborhood of the 
frequency 0=1. Also the dispersion curves for the 
modified plasma mode and the modified extraordinary 
electromagnetic mode approach each other again 
in the neighborhood of the frequency 0=7?. It may 
therefore be anticipated that the dispersion in the 
axial boundary layer might undergo rapid changes 
in the neighborhood of the two frequencies 0=7? 
and 0=1. 

3. Dispersion in the Axial Boundary Layer 

For the purpose of obtaining the proper approxi- 
mations in the axial boundary layer, it is convenient 
to start by rewriting (12a-d) as follows: 

A 10 =-C 2 ^l 2 l'2(n 2 -l)(i} 2 -l-R 2 ) + R 2 (l 2 -l) 

+ ^{(0 2 -l) 2 -OW}] (43a) 

^20=6 Y o 4 O 2 [(O 2 -l)(O 2 -7i >2 ) + 7? 2 (/ 2 -l) 

+ 20 2 -g{0 2 -l-Z 2 7i> 2 }] (43b) 

^30=-6> 2 O 4 [O 2 -7? 2 +7? 2 (l-/ 2 )] (43c) 

and 

A =(0 2 -l)[0 4 -0 2 (2+7? 2 ) + l] (43d) 

The axial boundary layer may be easily shown to 
be defined by those values of 6 for which 



l 1 t ^ 2 



(44) 



where t^>0 and is of the order unity. First, it is de- 
sired to restrict attention to the neighborhood of the 
frequency Sl=R. For this purpose, let 



Q*=R*+a 



CI 



(45) 



where s is of the order unity. The substitution of 
(44) and (45) in (43a-d) and the rearrangement of 
the terms in powers of a 2 /Cl, yields the following 
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A — P 2 \ A(o)\Aii)2l L/|( 2 > —4-/4 (3) - L/IW — 

(46a; 

/l __ / n n i \ A (0) I a (i) _^ _L A (2) — 4- /I (3) ^ 

^120 — ' <)'< ^20T^2fl ^72 I -^20 p± r^-20 ^6 

/I — — *° 2 rt4 /I (0) 4_ 4 (i) -?4 _1_ A (2) ^L 



(46b) 
(46c) 



and 



where 



A>=4T +A±» %+AP %+A™ % (46d) 



Alf=2R 2 (l-R 2 ) 



R 2 S 



A$=R*(s-t-2) 

A^=RH+R*s 

Ai a) = -2R i +3R 2 -l. 



(47a) 
(47b) 
(47c) 
(47,1) 



The expressions for all I lie other A's arc not needed in 
this analysis and hence are not reproduced here. 
The exact dispersion equation (11) maybe rewritten 
with the help of (46) and (47) as follows: 



[Ar+A , , ,, -S+^ 2, 7^ 



o C aJX 

P 2 A (0) 4- A < 1} - \-AW--\-A (3) — 4- /4 (4) — \\ 

L ( o ' o < o ' oj 



r„ 



(48) 



Note that all the A's are of the order unity. 

If the phase velocities of the order a or lower are 
only sought, it is legitimate to retain on the left side 
of (48) only the leading terms which are of the order 
Cla* and neglect all the other terms which are at least 
of the order a 2 jCl lower than the leading terms, 
resulting in the following quadratic equation in X 2 : 



A (0) A (0) 

a a A (0) a -\-a A (0) 







(49) 



provided ^L^^O. Using (47) in (49), it is found 
that 



(«)'+ 2 "(«)' +h '=° 

where 

«(R 2 -\)-R 2 (t + 2) 



(50 a) 



U 



4(/? 2 -l) 



n 2(R 2 -1) { } 



For t=0, which corresponds to the case of propaga- 
tion along the static magnetic field, the two modes 
whose phase velocities are specified by (50a) are de- 
coupled as evidenced by the fact that (50a) factors 

into two linear equations in ( - J yielding the follow- 
ing solutions: 



ex- 
ex 



s 
~2 

R 2 



R 2 -\ 



for t=0 



for £=0. 



For t7*0j the solutions of (50a) are given by 
QY =-U±^U 2 ~-W for^O. 



(51a) 



(51b) 



(52) 



It may be easily shown that U 2 —W can never vanish, 
with the result, the two phase velocities specified by 
(52) are always different. 

The phase velocity given by (51a, b) and by (52) 
for t = 1 are sketched in figures 9a and 9b for (i) R 2 <^1 

and (ii) A >2 >1, respectively. For R 2 <^lA-\ does 
not lead to real phase velocities and ( ) corresponds 

to real phase velocities for s<^—RH. Also ( - ■ Wo 

for s=-RH. For R 2 >1, (51a) for ,s<0 and (51b) 
for all s leads to real phase velocities. Also (X/a)i 
in (52) for all s and (X/a) 2 for s<l — R 2 t leads to real 
phase velocities. Also (\/fl)i>(\/d) 2 for s<^—RH 
and (X/a) 2 = for s= -RH. From (44) and (45), 
it is clear that s=—RH corresponds to tl=Rl. For 
R 2 ^>1 and £ = 0, the two phase velocities are found 
to be equal for s = —2R 2 jR 2 —l. 

In the axial boundary layer, in order to understand 
the physical significance of the variation with respect 
to s of the phase velocities of the order of t he acoustic 
velocity a in the electron gas or lower, for the fre- 
quencies in the close neighborhood of il = R, it is 
necessary to compare figure 9a with that portion of 
figures 2b and 8a enclosed by a solid line. The 
portion of the dispersion curve enclosed by a solid 
line in figure 8a is just the modified extraordinary 
electromagnetic mode, MEM (A! 11 ), whose phase 
velocity, according to (34), goes to zero for Sl=Rl. 
The phase velocity in figure 9a, according to the fore- 
going discussion, goes to zero for s=—RH, or £l=Rl. 
Therefore the phase velocity diagram for t^0 in 
figure 9a is just the expanded version of the dispersion 
curve of the MEM(A" n ) mode near Sl—R and \=a. 
When£ approaches zero, figure 9a shows that the dis- 
persion curve for the MEM(Z n ) mode just shifts to 
the higher frequency side and has its zero now at 
s=0 or Q,=R. A comparison with figure 2b shows 
that the phase velocity diagram in figure 9a for the 
for the case £=0, is the same as the purely transverse 
extraordinary electromagnetic mode propagating 
along the direction of the static magnetic field. 
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-20 



1/2 



Figure 9a. Dispersion in the axial boundary layer in the neighborhood of fi=R and X=a (R 2 =l/2). 



Summarizing the above results, it follows that for the 
case R 2 <1, the modified extraordinary electro- 
magnetic [MEM(J n )j mode in the neighborhood of 
&=R continuously goes over to the uncoupled, 
purely transverse electromagnetic mode as the direc- 
tion of propagation approaches and finally coincides 
with that of the static magnetic field. 

The portion enclosed by the solid line in figure 2c 
has two dispersion curves. The phase velocity curve 
which goes to zero for tt=R belongs to the purely 
transverse extraordinary electromagnetic mode and 
the other curve which is nearly horizontal belongs 
to the purely longitudinal plasma mode. These 
curves refer to propagation in the direction of the 
static magnetic field. From (16), the phase velocity 
of the plasma mode for Q=R is seen to be equal to 

aR 
/ p 2 _ ' There are two intersecting curves for the 

phase velocity in figure 9b for the case £=0, which 
refers to propagation along the static magnetic field. 



The curve for \, obviously, is a part of the purely 
transverse extraordinary electromagnetic mode since 
it has a value zero for s=0 or tt=R. The other 
curve, namely that for \ 2 , as seen from (51a), has the 
same phase velocity as that of the plasma mode for 
tt= R and hence is a part of the dispersion curve for 
the purely longitudinal plasma mode. Note that 
these two dispersion curves intersect and hence have 
the same phase velocities for s= — 2R 2 /R 2 —l. The 
portion enclosed by the solid line in figure 9b has 
two, nonintersecting dispersion curves. Obviously, 
the one on the low frequency side is a part of the 
modified extraordinary electromagnetic mode, since 
it has zero value for tt=RL The other curve, which 
lies on the higher frequency side, is a part of the 
modified plasma mode. These two nonintersecting 
dispersion curves refer to propagation in a direction 
which makes a small angle with that of the static 
magnetic field. On comparing the two dispersion 
curves enclosed by the solid line with the two curves 
in figure 9b for the case t^O, it is evident that the 
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Figure 9b. Dispersion in the axial boundary layer in the neighborhood of fi = R and X=a (R 2 = 4/3). 



phase velocity curves for Xi and X 2 belong to the 
modified plasma mode and the modified extraordi- 
nary electromagnetic mode [MEM(X n )] respec- 
tively. ^ The results depicted in figure 9b show the 
transition that takes place in the axial boundary 
layer for the case i? 2 >l. From an examination of 
figure 9b, it is clear that as the direction of propaga- 
tion approaches and finally reaches that of the static 
magnetic field, the higher frequency side of the dis- 
persion curve for the modified plasma mode merges 
with the lower frequency side of the dispersion curve 
for the MEM(X n ) mode to form a continuous curve 
which belongs to the uncoupled, purely longitudinal 
plasma mode. On the other hand, the higher fre- 
quency side of the dispersion curve of the MEM (X 11 ) 
mode crosses that of the uncoupled plasma mode, 
merges with the lower frequency branch of the modi- 
fied plasma mode to form a continuous dispersion 
curve which obviously belongs to the uncoupled, 
purely transverse extraordinary electromagnetic 



mode. Simple evaluation of the generalized 
Poynting vector may be made to establish that the 
two modes whose dispersion curves cross each other 
for 2=0 near Q=R and \=a are uncoupled. Thus 
in the axial boundary layer, for the case i? 2 >l, a 
part of the decoupling of the plasma mode takes 
place in the close neighborhood of ti=R and \=a. 

It is now desired to examine the behavior of the 
dispersion in the neighborhood of the frequency 
= 1. For this purpose, let 



Q 2 =l+s- 



(53) 



where s is of the order unity. Employing the same 
procedure as before and making use of (44) and (53), 
it may be easily derived that the phase velocities of 
the order Op or higher are specified by the following 
cubic equation in X 2 : 
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(l + t+2s) 



(& 



s-B 2 (s+t+2) 
B 2 s 



(£)' 



(~i? 2 +l) 



=0. (54) 



For £ = 0, which corresponds to the case of propaga- 
tion along the static magnetic field, the three modes 
whose phase velocities are specified by (54), are 
decoupled as evidenced by the fact that (54) factors 
into three linear equations in (\/a) 2 yielding the 
following solutions : 



and 
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for t=0. 



(55a) 
(55b) 
(55c) 



For the case ^0, it may be verified that the discrimi- 
nant of the cubic equation (54) in (A/Co) 2 is negative 
and as a result all the three roots are real and unequal. 
The phase velocities given by (54) for t=l and 
those given by (55) are plotted in figures 10a and b 
for the case (i) R 2 <^1 and (ii) R*\>1, respectively. 
For R 2 <il, (55a) for s^>0 and (55b) yield real phase 
velocities whereas (55c) does not give rise to real 
phase velocities. Also, it may be proved that the 
number of positive roots of (54) is either one or two 
according as s is less than or greater than zero. 
Note that only positive real roots will yield real 
values for the phase velocity. For R 2 ^>1 (55a) for 
5^>0, (55b) and (55c) yield real phase velocities. 
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Figure 10a. Dispersion in the axial boundary layer in the 
neighborhood of 12=1 and A=Co (R 2 = l/2). 



Also the number of positive real roots of (54) giving 
rise to real phase velocities is either two or three 
depending on whether s is negative or positive. 

For R 2 <C1, a comparison of figure 10a with those 
portions of figures 2b and 8a enclosed by a broken 
line shows that for the case £=0, which corresponds 
to propagation parallel to the static magnetic field, 
the curves marked 1 and 2 in figure 10a are the 
portions of the dispersion curves for the purely 
longitudinal plasma (P) mode and the entirely 
transverse ordinary electromagnetic (EM) mode 
respectively. The phase velocities of these two 
modes are equal for s=R/R-\-l. Similarly for t==l 9 
which corresponds to a direction of propagation very 
slightly inclined to that of the static magnetic field, 
the curves labeled 3 and 4 in figure 10a are the sec- 
tions of the dispersion curves belonging to the 
modified ordinary electromagnetic mode and the 
modified plasma mode respectively. It follows then 
from figure 10a that as the direction of propagation 
approaches and finally coincides with that of the 
static magnetic field, the branch of the modified 
ordinary electromagnetic mode for s^>i?/i?+l joins 
with that of the modified plasma mode for s<^R/R-\- 1 
to form a continuous curve which belongs to the 
uncoupled, purely transverse, ordinary electro- 
magnetic mode and the branch of the modified P 
mode for s^>R/R+l joins with that of the modified 
ordinary EM mode for s^R/R+1 to form a smooth 
curve which is part of the dispersion curve of the 
uncoupled, purely longitudinal plasma mode. It 
maybe easily verified that \/C =l is a solution of (54) 
for s= 1 for any value of R 2 or t. From this fact, it is 
evident that the modified plasma mode has its phase 
velocity equal to the free space velocity of electro- 
magnetic waves for the frequency given by £2 2 =1 + 
a 2 /Cl 
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Figure 10b. Dispersion in the axial boundary layer in the 
neighborhood of 12=1 and X=C (R 2 = 4/3). 
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For the case /? 2 >1, a comparison of figure 10b 
with that portion of figures 2c and 8b enclosed by a 
broken line enables the understanding of mechanism 
of coupling of the various modes in the axial bound- 
ary layer in the close neighborhood of 12=1. From 
figure 10b it is seen that as the propagation direction 
is changed from that of the static magnetic field, 
the branch of the dispersion cnrve for the plasma 
mode for s<^R/R-\-l joins with that of the ordinary 
EM mode for s^R/R + l to form a smooth curve 
which belongs to the modified ordinary EM mode. 
And the branch of the dispersion curve of the ordi- 
nary EM mode for s^R/R+1 joins with that of the 
plasma mode for s^>R/R+l and which in turn 
merges with that portion of the dispersion curve of 
the extraordinary EM mode for s^>R/R—l to form 
the dispersion curve corresponding to the modified 
plasma mode. In a similar manner the branch of 
the dispersion curve of the extraordinary EM mode 
for s<CR/R—l merges with that of the plasma mode 
for s^>R/R—l to form a smooth dispersion curve 
belonging to the modified extraordinary EM mode. 
From the foregoing discussion, the manner in which 
the plasma mode gets decoupled in the axial bound- 
ary layer in the close neighborhood of 12=1 and 
X=(7o, i s niade clear. 

With the understanding of the behavior of the 
dispersion of the three modes in the axial boundary 
layer, it is evident that the whistler mode corre- 
sponds to the modified extraordinary electromagnetic 
mode for all directions of propagation and no incon- 
sistency in the nomenclature of the whistler mode 
arises as in the magneto-ionic approximation. 

4. Concluding Remarks 

It is appropriate to make the following remarks in 
conclusion. The diagrammatic representation in the 
12 2 — R 2 space used in this paper for elucidating the 
regions of propagation of the various modes is very 
similar to the one first introduced by Clemmow and 
Mullaly [1955] and extensively used by Allis, Buchs- 
baum, and Bers [1962]. The two parameters used 
by these authors both contain the source frequency. 
Only one of the two parameters used in the repre- 
sentation of the parameter space employed in this 
paper, contains the source frequency and this feature 
facilitates greatly the determination of the propaga- 
tion characteristics explicitly as a function of 
frequency. 

The simple explanation that the compressibility of 
the plasma merely introduces a lower cutoff in the 
phase velocity curves obtained as a result of the 
magneto-ionic approximation, is frequently advanced. 
From the results of this paper, the magneto-ionic 
theory is seen to yield zero phase velocity for the three 
frequencies, 12=0, 12 3 and 12 4 . The compressibility 
of the plasma does not affect the phase velocity 
near 12=0 and merely shifts the zero near 12 3 to the 
high frequency side with the result that the zero 
phase velocity now occurs at 12=7*7. Only the phase 
velocity near 12=12 4 , instead of going to zero has a 



cutoff value equal to the acoustical velocity in the 
electron gas. Consequently the frequently advanced 
simple explantion of the effect of compressibility is 
not fully correct. 

Kieburtz [1963] has investigated the same problem 
that is treated in this paper using the coupled-mode 
theory and has found, for the range of frequencies 
given by 0<ci?<O in which the whistler phe- 
nomenon is observed, that there is coupling between 
the extraordinary electromagnetic mode (whistler 
mode) and the purely longitudinal plasma wave. 
The results obtained in this paper show that for 
R<Cl, there is no coupling between the whistler mode 
and the purely longitudinal plasma wave. There- 
fore, the results reported by Kieburtz appear to be 
incorrect. 

The frequency 12 4 , in the close neighborhood of 
which there is a sharp transition of the character of 
the wave from being predominantly transverse to 
being predominantly longitudinal is referred to in 
the literature [Allis, Biichsbaum, and Bers, 1962] as 
the frequency of plasma resonance. The results of 
this investigation show that for i?<0, the frequency 
of the plasma resonance is equal to the plasma fre- 
quency co„ in the direction of the static magnetic 
field and increases gradually to the value Va^ + co;; 
for the propagation across the static magnetic field. 
On the other hand, for 7?>1, the frequency of the 
plasma resonance is again equal to the plasma fre- 
quency cc p for propagation along the static magnetic 
field, "jumps" to the value of the gyromagnetic 
frequency co c as soon as the propagation direction is 
changed ever so little from that of the static mag- 
netic field and increases gradually and attains the 
same value V<4 + W c for propagation across the static 
magnetic field. 

An important and a satisfying result of this in- 
vestigation is that the magneto-ionic theory is sur- 
prisingly accurate over a wide range of t he parameters 
of interest. Where it fails and where inconsistencies 
set in as a result of its applicat ion, t he present theory 
may be used to obviate the difficulties. 
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